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In the recent years, field theory on non-commutative (NC) spaces has 
attracted a lot of attention. Most literature on this subject deals with per- 
turbation theory, although the latter runs into grave problems beyond one 
loop. Here we present results from a fully non-perturbative approach. In 
particular, we performed numerical simulations of the Xcj)* model with two 
NC spatial coordinates, and a commutative Euclidean time. This theory is 
lattice discretized and then mapped onto a matrix model. The simulation 
results reveal a phase diagram with various types of ordered phases. We 
discuss the suitable order parameters, as well as the spatial and temporal 
correlators. The dispersion relation clearly shows a trend towards the ex- 
pected IR singularity. Its parameterization provides the tool to extract the 
continuum limit. 
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1. Field Theory on a Non-Commutative Space 

The simplest way to introduce non-commutative (NC) coordinates is to 
impose the relation 

[x^,Xu] = iQ^u , (1-1) 

O being a constant, anti-symmetric tensor, while are Hermitian coordi- 
nate operators (which cannot be diagonalized simultaneously). This relation 

* Talk presented by W.B. at Workshop on Random Geometry, Krakow, 2003 
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leads to an uncertainty in c-space, Ax^Axj^ > {fi ^ u). The pre-history 
of this idea involves private communication among celebrities like Heisen- 
berg, Peierls, Pauli and Oppenheimer. The latter asked his student Snyder 
to work it out, which yielded the first publication about physics on a NC 
space followed immediately by Ref. ||2j. 

The mathematical framework for the formulation of field theory on NC 
spaces was worked out in the eighties However, a real boom of interest 
was triggered only in the late nineties by the observation that string and M 
theory at low energy in a magnetic background field can be identified with 
NC field theory JT. ^ This boom persists up to now, as is manifest from 
new preprints on this subjects appearing day after day. 

The identification of Refs. 5, 4 transforms the background field into 
a term. The same idea has been around in solid state physics long be- 
fore, where it led to a new description of the quantum Hall effect [H]. The 
origin of this idea goes back to Peierls: an electron in a strong magnetic 
background field has an obvious description by NC coordinates, where the 
magnetic field is replaced by non-commutativity, with an extent inverse to 
the magnetic field strength. 

To illustrate this fundamental property in the simplest possible way, we consider 
an electron moving in a plane, with position x = {xi, a;2, 0), exposed to a magnetic 
field B = (0, 0, B). If this field is strong, the Lagrangian 

L[x,x\ ^ —x'^ + eBeijXiXj , (i,j = l,2) (1.2) 

can be reduced to the second term, so that the canonical momentum reads 



dL 

IIj = ^T- = eBeijXi . (1.3) 



Applying now the canonical quantization rule 

[xt,Uj] ^ ihS.j (1.4) 

we arrive at 

[x^^Ili] = eBeij[xi,Xj] . (1.5) 
Indeed, together with cq. H1.3I) this corresponds to the NC relation 

[x,,Xj] = iOij , Qij ^ -^^Cy ■ (1-6) 

This illustration of Peierls' map, along with a more precise derivation based on the 
Hamiltonian formalism, is explained for instance in Ref. jTj- 



^ Strictly speaking, also that observation occurred in the literature much earlier 
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Hence one motivation to study NC field theory is simply its application 
as a formalism to describe certain effects in the commutative world. Such 
effects are typically related to a background field, which is then transformed 
away by going non-commutative. 

However, in addition to that concept, the present fashion also includes 
studying the possibility of a really NC space-time. A deep, qualitative 
difference from the commutative space-time is the occurrence of a non- 
locality of range \/||0j[. Obviously this feature raises conceptual problems, 
but from the optimistic point of view it is a source of hope for a link to 
quantum gravity. 

In fact, there is a claim that attempts to merge quantum theory and 
gravitation imply quite generally a NC space. To illustrate this line of 
thought, we quote a simple Gedankenexperiment. 

Some event is measured with accuracy Ax, Ay, A^, At. This requires an 
energy concentration, which implies a gravitational field. In the extreme case, 
the latter imposes an event horizon beyond the uncertainty, so that the event is 
effectively invisible. One may now evaluate the condition for avoiding this, i.e. for 
dealing with actually detectable events. On the Planck scale, an estimate in Ref. 
[H] suggests the constraints 

AxAy + AxAz + AyAz > 1 , 

{Ax + Ay + Az)At>l, (1.7) 

so that the NC space seems indeed natural as soon as gravity is involved. 

However, this argument should come along with at least one remark of cau- 
tion: much of the literature excludes time from the non-commutativity, including 
the remainder of this article. Otherwise the problems related to causality 1^ are es- 
pecially severe. From the above argument, however, this step would not be justified. 

The last point is also related to the question if and how the Wick rota- 
tion from an Euclidean to a Minkowski signature can be performed in the 
NC world. This is not ultimately settled yet, but in our case of a commuta- 
tive time the issue is less problematic. Anyhow, here we work in Euclidean 
space, and we are happy there, without worrying about the details of the 
transition to the NC Minkowski space. 

Since the idea that our space could really be NC is fashionable, of course 
there are already numerous speculations about possible measurements of Q. 

One suggestion is based on the deformation of the photon dispersion relation 
due to 6 |10j. Blazers (highly active galactic nuclei) emit bursts of photons over a 
broad energy spectrum. Assuming this emission to be simultaneous, a relative de- 
lay could in principle establish bounds on O. However, in addition to experimental 
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difficulties, the knowledge about the deformed dispersion is also limited to a one 
loop calculation IIJJ. 

The last limitation is especially worrisome in the light of the fact that 
most higher loop calculations are not feasible yet — no systematic ma- 
chine is known for them. Perturbation theory is even more complicated 
than in the good old commutative space, which is in striking contradiction 
to the original hope that G would simplify the perturbative renormaliza- 
tion It is true that part of the UV singularities are removed due to 
the non-commutativity, but others remain, in particular those in the planar 
diagrams • What makes the situation worse is that the non-planar diver- 
gences do not just disappear, but they are rather turned into IR singularities 
with respect to external momenta. This effect is denoted as UV/IR mix- 
ing |13j . At this point we want to give again just a simple intuitive reason; 
we will be somewhat more explicit below in the framework of the Xcp^ model. 

We return to simple natural units h= 1, without involving the Planck scale any 
more. If we combine Heisenberg's uncertainty Axj ^ ^/^Pj with the NC relation 



we see that for Apj — > the other momentum components explode, Api — > oo, and 
vice versa. This also suggests that in addition to the Heisenberg term, Axj picks 
up a term linear in the Api, which may be denoted as a "string modification" of 
the uncertainty principle. 

In the work to be presented here we are going to consider a 3d Euclidean 
space with a commutative Euclidean time t and two NC spatial coordinates, 
which obey 



We are happy to avoid the horror of NC perturbation theory by taking the 
fully non-perturbative approach of numerical simulations. To this end, the 
space should first be lattice discretized, as in the commutative world. Here 
we also need a second step, namely a mapping of the lattice field theory 
onto a matrix model, which will be described in Section 2. At this point 
we only comment on the general structure of a 2d NC lattice of spacing a, 
following Ref. [Hj . 

The restriction of the spectrum of to the lattice sites corresponds to the 
operator identity 



Ax J Qij/Axi QijAp, 



^ j) 



(1.8) 



[x^,x^] = i6e^ 



(1-9) 




(1.10) 
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As in the commutative case, we want the momentum components fc^ to be periodic 
over the Brillouin zone, 

exp(i/c^£p) . (1-11) 





27r- 






a . 





exp yi 

Multiplying both sides by the factor exp(— ifc^a;^) now leads to consistency with 
eq. iff 

G ^ - (1-12) 

which has amazing consequences: any NC lattice is automatically periodic, say 
over a lattice volume L x L. Then the discrete momenta fc^") = occur, where 
rip G Z. The non-commutativity parameter can be identified as 



e^-a^L. (1.13) 

TT 

We see that the continuum limit a — s- and the thermodynamic limit L oo 
are manifestly entangled, which is again an aspect of UV/IR mixing. Taking these 
two limits simultaneously in such a way that 9 remains constant is denoted as the 
double scaling limit. 

2. The non-commutative Xcf)'^ model 

NC field theories can be formulated in a form which looks similar to the 
commutative world, if all the fields are multiplied by the star product (or 
Moyal product), 

(\ d d \ 
f{x)-kg{x) •= (^2^®^''9^9^ ) f^^'^3{y)\x=y ■ (2.1) 

In the particular case of bilinear terms in an action, the star product is 
equivalent to the ordinary product, hence in these terms the star product 
is not needed. 

Based on these rules, we can write down for instance the action of the 
NC Xcj)^ model, 

S[ci>\ = j d''x\^d^(t>d^<P+'^4>'' + ^cp^c^^c^*c^^ . (2.2) 

Since only the self- interaction term involves 0, the coupling strength A also 
determines the extent of NC effects in this model. 

To render the above star product rules plausible, we consider the composition 
of plane wave operators, 

g.p.s. . g^<?„x„ ^ _ '-p^,Q^,uq.) e'(f+«)-*- . (2.3) 
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This provides a prescription for the translation from the commutative product 
gipx^iqx ^^(j formulation in terms of x, i.e. without using the operator x, 

^ipx^iqx _^ ^-■^P^&^.^qi. gi(p+q)x ._ ^ipx ^ ^iqx ^2 4) 

More generally, we can adopt this translation rule to the product of fields, which 
can be decomposed into plane waves, 

(t){x)il}{x) (t){x) ei^A'^A'-^- ^]{x) := 4>{x) * V'(x) . (2.5) 
Turning now to the peculiarity of bilinear terms, it is easy to see that in 

(fx(l){x)i^ip{x)^ J d'^xct){x) 1~ ^df,ef,^dl + ... ^p{x)^ J d'^x(j}{x)^ix) 

all the terms in the square bracket which involve Q cannot contribute, based on 
partial integration and Q^^ — —Qu^. 

The one loop level of this model is suitable for the illustration of the 
general properties of NC perturbation theory that we mentioned in Section 
1. 

To this end, we consider the one loop level of the one particle irreducible 2-point 
function to the action 1)2. 2|l . i.e. the n = 1 contribution to 

CXD 

-^A"r(")(p) . (2.6) 

It contains a planar and a non-planar term, 

^(1) - i /" '^'^fe ^ p(i) _ 1 /" '^'^^ exp{ikf,Qf,i,p„) 



planar 3 J (^^Tr)^ + ' non-planar g J (^^Tr)^ + 

(2.7) 

which are illustrated in Fig. We see that the planar part is independent of Q 





planar non-planar 

Fig. 1. The planar and non-planar one loop contribution to the two-point function 
^2.b]) . given in eqs. \2.'l]) . 

|12| . whereas the non-planar part is affected by UV/IR mixing ^3]. To reveal what 
this means, we introduce a momentum cut-off A and obtain I14| ( for m > 0) 



(1) ,^ m(''-2)/2 4 N(2-<i)/4 . n ' ^ 

^non-planariP) = g (2^)rf/2 ( ~ Pl^^t^.P'^ } [m^ — - p^Bf^.p, ^ 
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where = Q^pQpi, and K is the modified Bessel function. In particular in d = 4 
the divergent part is given by 




A5- — PtJ.'^p.uPy 



In general the effective cut-off Acff remains finite in the UV limit A — > c», but it 
diverges if we take simultaneously the IR limit p — > 0. 

A comprehensive one loop study of the NC Xcf)^ model in various dimen- 
sions has been performed in Ref. That work dealt with a self-consistent 
Hartree-Fock type approximation, which would be exact for the 0{N) model 
at large A^. The authors assumed its validity also for the A^ = 1 theory and 
conjectured in particular a prediction for the phase diagram. As a general 
feature, the system undergoes some ordering at m? <C 0, which corresponds 
in some sense to a very low temperature. If rm? is lowered to that point, 
Gubser and Sondhi predict in d = 3 and d = 4 the following behavior: 

• At small 9, there is an Ising transition to a uniform order, as in the 
commutative case. 

• At larger 6, the ordered state has a structure of stripes or even more 
complicated patterns. 

Note that the formation of such a stripe order implies the spontaneous 
breaking of translation invariance. Therefore Gubser and Sondhi did not 
expect this phase in d = 2. We will return to this point at the end of Section 



The same question was also studied by means of a renormalization group 
analysis [TB] . 

Our goal was a non-perturbative, quantitative verification of that quali- 
tative conjecture, and our methods are Monte Carlo simulations. The lattice 
formulation (c.f. Section 1) is not hard to write down, but it is very hard 
to simulate due to the star product. The way out which enabled efficient 
simulations is a mapping onto a dimensionally reduced matrix model, as 
suggested in Ref. ilX,. This method is a refinement of the matrix model 
approach to NC gauge theory in the continuum jl8j . 

Assume the field cl){x, t) to be defined on a A^^ x T lattice of unit spacing. 



3. 
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According to Ref. jl7j the following action is equivalent to the lattice action: 

T 2 



t=i ^=1 

+^(^4>{t+i)-^{t)y+^p{t)+^4>Ht)] . (2.9) 

Here (^(t) (t = 1 . . . T) are Hermitian N xN matrices. In the (commutative) 
time direction the kinetic term takes the ordinary discrete form. The quar- 
tic term — which is in the original formulation plagued by repeated star 
products — looks relatively simple now, although it should be noted that 
these matrices are not sparse (in contrast to typical lattice formulations of 
commutative field theory), hence the fourth power requires some computa- 
tion. However, much of the complication due to the NC geometry is now 
manifest in the non-standard kinetic term in the spatial directions. The 
so-called twist eaters provide the shift which corresponds to one lattice 
unit, if they obey the 't Hooft algebra 

r^^r^ = z^^r^r^ . (2.10) 

The tensor Z^^, = Z*^ is called the twist factor. In general it may have the 
form Z12 = e^'^**^/^, where k ^Wj. The representation of the twist eaters 
that we choose |2U1 121j requires the twist to have the specific form 

Z12 = exp [^^^^) , (2.11) 

as we are going to sketch below. Comparison with the general form shows 
that then N has to be odd. The lattice model and the matrix model are 
connected by Morita equivalence, which means that their algebras are fully 
identical. 

Let us go back to the NC lattice formulation discussed in Section 1. Since the 
momenta are discrete, we should not use the (unbounded) operators x^, /i = 1,2 
to describe the lattice sites. Instead we introduce the unitary operators = 
exp{j^ix f^) , which obey the commutation relation 



Z^Z^ = exp ( - ■^^iSf.S^Zi.Zfj, = exp - ^ie^i/ ) i'lyi'p , (2-12) 



where we used = and eq. (|1.13() . The lattice translation operator Z)^ = 
exp(a9;^) is supposed to fulfill the relation 

Df^Z^Dl = e^^'^^^-Z^Z^ . (2.13) 

^ Its standard form, which occurs for instance in the "twisted Eguchi-Kawai Model" 
|19|. uses k — 1. 
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The issue is now to find a matrix solution for the conditions H2.12|l and H2.13(l . This 
solution is unique only up to symmetry transformations, hence one may end up 
with different twist factors. We assumed the unitary twist eaters to take the form 



Ti = 



/ 





V 1 



1 



1 

/ 



r. = 



/ 1 










221 







Z|i . 






^21 




V 




■ ) 



where Z12 — Z21 characterizes the twist. This ansatz solves indeed the conditions 
(|2.12(l and H2.13|l . iff we choose the twist of eq. H2.11(l 21 . However, a different 
may lead to an alternative solution with a modified twist, 



ansatz for D 
see e.g 



^ resp. 
Ref. 122 ■ 



3. Numerical results 

An overview of our numerical results is given by the phase diagram in 
Fig. [21 For a suitable re-scaling of the axis, this diagram stabilizes with 
increasing N = T. At strongly negative m? (corresponding to a very low 
temperature) the system seeks some order. At small self- interaction A the 
non-commutativity is not much amplified and the systems behaves as in the 
commutative world, i.e. it is ordered uniformly. At larger A a new striped 
phase sets in, in the spirit of Ref. J^l- Such a phase is not known in the com- 
mutative A(/>^ model, but similar phenomena appear in solid state physics, 
see for instance Ref. [23j. In that phase, a non-zero mode condenses, so 
that the ground states correspond to a stripe pattern, or to more compli- 
cated checker- field- type patterns. Typical examples for the various cases 
are shown in Fig. |21 

The transition from disorder to order can be localized well and it appears 
to be second order. On the other hand, the transition between the uniform 
and the striped phase is more difficult to localize and it is expected to be 
of first order. 

The order parameter that was used here is based on the spatial Fourier 
transform of the scalar field, (f){p, t). We average over the time direction and 
turn it such that an eventually condensated mode can be detected optimally. 
This is done be introducing the function 



M{k) - ^ "^^^ 



J^rp k=\p\N/2TT 



t=i 



(3.1) 



and its expectation value is our order parameter, sensitive to the mode k. 
In particular, (M(0)) is the standard order parameter for the Z2 symmetry, 
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eg 







100 200 300 400 500 600 700 800 



Fig. 2. The phase diagram of the NC model, explored by means of numerical 
simulations. Since it contains a new type of phase we conclude: this diagram of 
phases simply amazes. 

(M(l)) is the staggered order parameter that detects patterns of two stripes 
parallel to one of the axes as in Fig. |31below on the left, (M(\/2)) detects a 
diagonal two-stripe pattern as in Fig.|31below in the center, (M(2)) detects 
a pattern of four stripes parallel to one of the axes as in Fig. |31below on the 
right, etc. We also measured the connected two-point functions of M{k), 
the peak of which allowed us to localize the order-disorder phase transition 
in the diagram of Fig. I^Jto a high precision. 

Next we consider the spatial correlator at a fixed time. 



As it can be expected, in the disordered phase this decay is fast, but in 
contrast to our daily experience in the commutative world, it is in general 
not exponential (nor polynomial), as Fig. 0] (on the left) shows. It may be 
plausible that the spatial decay behavior is somehow distorted by the NC 
geometry; what is perhaps more surprising is that the decay moves closer 
to an exponential as A increases (so that the NC effects are expected to be 
enhanced), see Fig. |1] (on the right). This observation might be explained 
in the light of the pole structure described in Ref. [T^ . 

Turning now to the striped phase, the pattern with stripes parallel to 
one of the axes can be visualized well by plotting C(xi,0) and C{0,X2). In 
Fig. [21 we show two examples where the cuts through a configuration of two 
stripes parallel to an axis (on the left) resp. two diagonal stripes (on the 
right) can easily be recognized. 

Now we want to consider the correlations in the direction of the Eu- 
clidean time. More precisely, we consider the correlators of two fields (j){p, t) 




(3.2) 
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Fig. 3. Snapshots of typical, well thermalized configurations in different sectors of 
the phase diagram: disordered and uniformly ordered (above), and patterns of two 
stripes parallel to an axis, two diagonal stripes and finally four stripes (below). The 
dark and bright areas correspond to (j) > 0, (f) < 0. 

at the same momentum phut at different times. For instance, Fig. |Hlon the 
left shows 

G(t) :={4>{p = 0,t)^{p = 6,t + T)) (3.3) 

in the disordered phase, close to the phase transition. Since the time direc- 
tion is commutative, we rather expect here the usual behavior, and indeed 
these correlators follow neatly the exponential with periodic boundary con- 
ditions, i.e. a cosh shape. This shape allows us to extract the energy at rest 
— or effective mass — in the exponentially dominated sector (not too close 
to the center) as 

E{p = 6) = -ln^^^-^ . (3.4) 
G(t) 

By varying r we find a convincing plateau, see Fig. El on the right. 

We can now repeat this procedure also at non-vanishing momenta p, 
which finally yields the full dispersion relation E{p). We still stay in the 
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(a) A = 0.06, = -0.015 



(b) A = 2, 771^ = -0.3 



Fig. 4. Decay of the spatial correlator C{x) defined in eq. 13. g|) at N = 35 in 
the disordered phase, close to the ordering phase transition. On the left we see a 
clear deviation from the exponential decay. At larger X (on the right) the standard 
exponential behavior is approximated. 
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(a) A = 0.6, = -0.7 



(b) A = 6, 



Fig. 5. The profile of C{xi, 0) and C(0, 0:2) at N = 35 for two stripes parallel to an 
axis (left) and two diagonal stripes, as they occur at stronger coupling (right). 



disordered phase, close to the ordering transition. If we are close to the 
uniform order (small A), the dispersion has its usual linear shape (up to 
finite size corrections) as in the commutative case, see Fig. |3 on the left. 
However, the shape changes if we move to the vicinity of the striped phase, 
see Fig. EJon the right. Now the rest energy increases, since the effects of 
UV/IR mixing set in. In the large N double scaling limit it is expected to 
diverge. The minimum moves to a finite value of the momentum, which 
corresponds to the stripe patterns shown in Fig. |21 If we enlarge the lattice, 
we expect the dispersion relation to stabilize if the axes are taken in physical 
units (an exception is of course the close vicinity of p = 0). The required 
re-scaling of the axes determines the physical lattice spacing a. A clear hint 
that such a stabilization can indeed be achieved is the preliminary Fig. IS] — 
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Fig. 6. On the left: the temporal correlator G{t), defined in eq. ij^.^^) . at N = 35, 
A = 0.29, m? — —0.11. On the right we show the rest energy extracted at different 
values of r according to eq. which has a clear plateau. 
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(a) N = 45, X = 0.035, 
-0.009 



(b) iV = 55, A = 55, = -15 



Fig. 7. The dispersion relation E'^lj)'^) in the disorder phase, close to the ordering 
transition. On the left: at weak coupling — above the uniform phase — it takes the 
standard shape known from the commutative world, including the lattice correction 
(dashed line). On the right: at strong coupling — above the striped phase — the 
rest energy increases and the energy minimum moves to finite momenta. The line 
is a four parameter fit to the one loop prediction 12 U YMj . Although we are still in 
the disordered phase, the system already feels the trend towards a stripe formation. 



for the corresponding definition of the a (which also stabihzes the spatial 
correlator) we refer to Ref. j2H|- We see here in particular that the energy 
minimum stabilizes in physical units, which means that at large N we find an 
infinite number of stripes with a stable average width, which can be denoted 
as a zehra pattern. (Of course this also includes the interference of stripes 
in several directions; this is what we meant by checker-field pattern.) Once 
this is fully demonstrated, we have established the continuum limit and the 
final confirmation of the existence of the striped phase, as conjectured by 
Gubser and Sondhi fTSj . 
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Fig. 8. The dispersion relations at different sizes N and physical lattice spacings a. 
The product Na^ — 100 is kept constant, which corresponds to the prescription of 
the double scaling limit. We see that the rest energy seems to diverge at large N , 
whereas the shape of the rest of the dispersion relation stabilizes. 



We finally remark that the occurrence of stripes in the lattice formulation was 
also observed in the case of two NC dimensions in d = 2 |22l 1201 12] and in d = 4 
|25| . In both cases, this observation is not trivial from the theoretical point of 
view: since such a phase breaks the translation invariance spontaneously, one may 
wonder if it can possibly exist in d = 2 at all ^S], c.f. Section 2. However, there 
is no contradiction with the Mermin- Wagner Theorem because the proof of the 
latter assumes locality and an IR regular behavior, which is both not provided 
here. In d = 4, on the other hand, the situation might be special because we are 
dealing with the critical dimension in view of the renormalization group |16| . and 
an exciting question would be if there are any news about triviality once O enters 
the game. 



4. Conclusions 

We have presented results from a numerical study of the A(/>^ model in 
3d Euclidean space, where the two spatial directions are NC, whereas the 
Euclidean time is commutative. The system was first lattice discretized and 
then mapped onto a dimensionally reduced matrix model. On each time lat- 
tice point one obtains a Hermitian N x N matrix, where N xN is the spatial 
lattice volume. We denote the physical (i.e. dimensionful) lattice space as 
o. The double scaling limit — > oo, at Na"^ = constant describes the 
limit of an infinite volume in the continuum at a finite non-commutativity 
parameter 9. This means that the UV limit and the thermodynamic limit 
are intertwined in this case, as a consequence of UV/IR mixing. 
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We find a phase diagram which stabilizes for > 25. The ordered regime 
sphts into a uniform and a striped phase, which is consistent with the con- 
jecture by Gubser and Sondhi. We discussed the order parameter, which is 
able to detect different types of stripe patterns. The spatial correlators in 
the disordered phase decay in some irregular way (fast but in general not 
exponential), whereas the correlators in time direction — at fixed spatial 
momentum — do decay exponentially. Based on the latter property we can 
determine the dispersion relation, which reveals the UV/IR mixing again. 
At large A the energy minimum drifts away from zero, in accordance with 
the occurrence of stripe patterns in the ground state. The 0-deformed dis- 
persion relation at large A seems to be IR divergent. Hence the short-ranged 
non-locality has once again also long range manifestations. 

As an outlook we hope to identify a similar dispersion relation for the 
photon, which would be of immediate phenomenological interest, as we out- 
lined in Section 1. We repeat that the experimental search for bounds on 6 
based on a possible frequency dependence of the speed of light rely so far 
on a one loop calculation TP (with an uncertain behavior at higher loops). 
Therefore a non-perturbative result would certainly be valuable. A first step 
in that direction was the study of Wilson loops in the 2d NC U{1) gauge 
theory 
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